Hydrodynamics of an ultra-relativistic fluid in the flat anisotropic cosmological model 
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Motion of an ultra-relativistic perfect fluid in space-time with the Kasner metrics is investigated 
by the Hamiltonian method. It is found that in the limit of small times a tendency takes place 
to formation of strong inhomogeneities in matter distribution. In the case of slow flows the effect 
of non-stationary anisotropy on dynamics of sound waves and behaviour of frozen-in vortices is 
considered. It is shown that hydrodynamics of slow vortices on the static homogeneous background 
is equivalent to the usual Eulerian incompressible hydrodynamics, but in the presence of an external 
non- stationary strain velocity fleld. 
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Dynamics of a relativistic fluid in a gravitational field 
has a set of distinctive properties in comparison with hy- 
drodynamics in a flat space-time [|| , ||] . An explicit de- 
pendence of the metric tensor components gik (t, r) on 
time coordinate t and spatial coordinates r = {x, y, z), to- 
gether with back influence of the matter on gravitational 
field, result in a wide variety in behaviour of solutions for 
equations of motion. Many phenomena in the general- 
relativistic hydrodynamics have not been studied yet. In 
order to get a better understanding of various processes 
taking place in the system, it is useful to analyse simple 
limit cases, when only one of many general-relativistic 
effects is mostly important. If a dissipation may be ne- 
glected, then one can simplify significantly the analysis 
by using the Hamiltonian formalism, application of which 
to hydrodynamics has given many interesting results (see 
the reviews Q and references therein). A brief dis- 
cussion of adaptation of the Hamiltonian method for the 
relativistic hydrodynamics can be found, for instance, in 
the papers and So, in the work j^, slow isentropic 
flows of an ideal relativistic fluid have been considered 
and a partial study of influence of a spatial inhomogene- 
ity on the motion of hydrodynamical vortices has been 
performed. On the contrary, in present work we suppose 
that the space is homogeneous and concentrate our at- 
tention on the effects of its non-stationary anisotropy. It 
is known that there exists the class of solutions of the 
Einstein equations for empty space-time, which has been 
found by Kasner in 1922, with metrics of the following 
form |]: 



V + * '(^-z), (1) 



where Ai , A2 , A3 are any three numbers satisfying two 
conditions 



Ai + A2 + A3 



0, 



Xf + Xj + Xl^ 2/3. 



(2) 



At small times, this metrics describes approximately also 
the case of Universe filled uniformly by a fluid, since in 
the limit t —^ a back reaction of matter on the gravi- 
tational field appears to be negligible ||^, |Q. In present 
work we have used this result and considered the hydro- 



dynamics of a perfect fluid having ultra-relativistic equa- 
tion of state, in the space-time with the fixed metrics (|^). 
We have established that at early times, fiows with a fi- 
nite mean momentum are in such a dynamic regime that 
a motion of each fluid element is nearly independent on 
the motion of other fluid elements. Due to this reason, 
there exists a tendency to formation of strong inhomo- 
geneities in the spatial distribution of the matter which 
is related to singularities of an approximate Lagrangian 
mapping. We have found also that even in the case of 
slow flows with zero mean momentum the nonstationary 
anisotropy influences strongly both on the dynamics of 
sound oscillations and on the vortex motion, making their 
properties be different in some aspects from properties of 
the analogous objects in the flat space-time. 

The plan of further exposition is the following. At 
first, a necessary explanation will be given concerning 
the employed theoretical method, and the expression for 
the Lagrangian C of relativistic ideal hydrodynamics in a 
given gravitational field will be written. After a brief dis- 
cussion of specific properties of corresponding equations 
of motion, we shall pass to the Hamiltonian description 
of the system dynamics and we will introduce the exact 
expression for the Hamiltonian Ti. corresponding to the 
ultra-relativistic equation of state of the matter. Some 
important properties of the dynamics in the limit i — > 
will be discussed. A quadratic (on small deviations) part 
Ti^^'' of the Hamiltonian will be used for the derivation 
of simplified equations of motion for sound waves and 
vortices. 

At first we introduce needed concepts. The energy 
density e, measured in the locally co-moving reference 
frame, and the density n of number of conserved parti- 
cles are connected by equation of state e = e{n) Q. The 
scalar n is equal to the length of the 4-vector of current 
— n{dx^/ds) Let us introduce the relative den- 
sity p{t, r) in such a manner that it obeys the standard 
continuity equation 



Pt+V{pv) = 0, 



(3) 



where 3-velocity field v(t,r) = = (77^,11^, w^) is de- 
fined as dx'^/dt on the world- line of the fluid point pass- 



1 



ing through {t,r). The equation (^) follows from the 
equation n*j = 0, expressing the conservation law for the 
fluid amount , if p and n are connected by the relation 



\/ 900 + '2goaV" + gaf3V"vf^, 



(4) 



where g = det\\gik\\ is the determinant of the metric ten- 



sor. 



The action functional S — J Cdt of relativistic hydro- 
dynamics is defined through the Lagrangian C = C{p, v} 
in the following way 



C = - 



-^goo + 2.9oaW" + gaf3v^vf>^ VWldr. (5) 



The equation of motion for the velocity field v{t, r) has 
the structure (generalized Euler's equation |^) 

This is merely the variational Euler-Lagrange equation 
expressing the least action S principle for variations of 
world-lines of fluid particles. The 3- vector p(t, r) = 
il/p)iSC/Sv), 



Pa 



e'{n) 



-{goa + ga/lV^) 



(7) 



is the canonical momentum of a liquid particle in the 
point {t,r). We sec that the definition of the canonical 
momentum in relativistic hydrodynamics depends in a 
complicated manner on the equation of state. The rela- 
tion (|^) can be resolved respectively to the velocity, which 
gives the inverse dependence: v = v{p, p}. The remark- 
able equation for the vorticity field n{r,t) = curlp(i, r) 



fit = curl[v X Jl], 



(8) 



can be obtained from the equation (|^) by applying the 
curl-operator. It implies that vorticity is frozen in the 
fluid and the generalized Kelvin's theorem (conservation 
of circulation of the field p along arbitrary closed contour 
j{t) transported by flow) holds: 



® (p • dl) = F-y = const. 



hit) 

Passing to representation of flows in terms of the fields 
p and p and defining the Hamiltonian Ti-ip, p} in such a 
manner: 



n{p,p} 



v=v{p,p} 

we have the equations of motion in the form 



(9) 



(10) 



Pt 



sn 



curl p 
P 



- V 



sn 



(11) 



In particular, for the potential flows, when the vorticity 
is equal to zero and p = Vyj, the dynamical variables p 
and Lp are canonically conjugated: 



Pt 



5n{p,^v} 

5if 



-ft 



5n{p,^f} 
Sp ■ 



(12) 



Equations (|l^) will be used further (in linearized form) 
for analysis of acoustic mode dynamics. 

Other interesting dynamic regime can take place if 
there is an equilibrium solution p = po(r), p = Po = 
const, independent on time. If vortical component of flow 
is small and the sound waves are excited weakly, then per- 
turbations of density are negligibly small in comparison 
with its equilibrium value. In such circumstances a slow 
dynamics of the vorticity is described by the equation , 



curl 



curl 



In 



n 



Po(r) 



(13) 



where the Hamiltonian of the vorticity 7i*{r2} is obtained 
from the quadratic (on Jp = p — po and Sp = p — po) 
part Ti^^^{Sp,Sp} of the total Hamiltonian after fixing 
the potential ip by condition of zero density variation 
Sn^^^lSip = 0: 



(14) 



It should be noted that although the equilibrium solu- 
tion must be independent on time, the Hamiltonian it- 
self may contain an explicit dependence on t. It is the 
case in the space-time with the Kasner metrics (^. In 
such a Universe the homogeneous equilibrium solutions 
p = Pq = const, p = Po = const are possible which 
fact can easily be seen from equations of motion (H) and 

Now we pass to concrete calculations. We will use the 
Kasner metrics (|l|) and the ultra-relativistic equation of 
state 



e{n) 



,4/3 



(15) 



since n — > oo when i ^ For convenience, we change 
definition of the exponents: 3Aq = 2/ic and perform the 
substitution t^/^ = r for the time coordinate in the action 
S, and also choose in an appropriate manner scales for 
spatial coordinates and for relative density p, in order to 
obtain 



(16) 



The procedure of calculation of the Hamiltonian 7i{p, p}, 
corresponding to this Lagrangian, includes first of all the 
solution of the relation p = {1/ p){S£/Sv), i.e. 
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(17) 



respectively to v, and its substitution to (^). Let us 
introduce the notations 



P 



We define the auxiliary function h{Q) by the conditions 
Q3(1-^)2^A3, h = ^AQ + ^ Af'"^ . (19) 

The exphcit expression h{Q) can be found with the help 
of formula for roots of cubical polynome. Designating for 
brevity 

1/3 

R{Q) = 2"i/3 (27 - 2Q3 + 3\/3v/27 - 4g3) , (20) 
we obtain that 



(21) 



We plot the dependence h versus y/Q in the Fig. |l|. 




FIG. 1. The dependence h on y/Q. 
At small Q, the approximate equality is valid 

Q < 1, 



and at Q ^ 1 the following behaviour takes place: 

1 

ig' 



KQ) 



> 1. 



(22) 



(23) 



As the result, the Hamiltonian of the ultra-relativistic 
fluid in the Kasner space-time has the form 



(24) 



For definiteness, we arrange the Kasner exponents in 
the order /ii < /j,2 < /J.3. The constants jia satisfy the 
inequalities 

- 2 < 2/ii < -1, -1 < 2^2 < 1, 1< 2/13 < 2, (25) 

following from the conditions (^. Let us suppose z- 
component of the mean momentum be different from 
zero: pz{) 7^ 0. We consider here the case when devia- 
tions of p from the mean value are comparatively small. 
Then at sufhciently small times the quantity Q reaches 
large values, since the order of magnitude of quantities 
p and p is conserved due to the conservation laws for 
fluid amount and for mean momentum. Asymptotics ( |2^ ) 
gives the approximate Hamiltonian in this limit 



pV4 



y ' 



•Pa 



jdr (26) 



If we neglect here the quadratic on p term, where the 
coefficient in front of p^ is small in the limit of small 
times, then equation of motion for the canonical momen- 
tum field does not depend on p at all. Each point of fluid, 
labeled by a label a = a", moves almost independently 
on other points in accordance with the approximate law 



x"(T,a) 



Pa{ai)dT 



Pa(a)ri+^3-2p„ 
|pz(a)|(l +P3- 2po 



(27) 



Deviations of arbitrary constants of motion p(a) from 
the mean value make possible finite-time singularities of 
the mapping r(T, a). Accordingly, quasi- two-dimensional 
domains of pancake-like shape may arise in the space 
where the density is much larger than the mean density: 
p = p(a)det||9a/9r|| po- There, the neglection of 
the quadratic on p term is no more valid, and actually a 
redistribution of momenta between fluid elements occurs. 

To investigate slow flows on a resting background, it is 
necessary to obtain the expression for p}. Let us 

expand Ti-lp, p} to the second order terms on Sp{T, r)_= 
/o(r, r) — Po and p(r, r). Then using of the expansion (j2^) 
gives 



2/3 

^(2) ^ Po_ 



--2Mc 



pldr 



Po 



2/3 



{Spf 



■dr. (28) 



Expression ( p8| ) and equations (|lj) allow us immediately 
to write the linearized equations of motion for the spatial 
Fourier- harmonics of sound waves: 



dr 



1 -1 
3P0 Vk, 



dpk 
dr 



PS 



-2^10 



fca^k. (29) 
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We can reduce this system to one ordinary linear differ- 
ential equation 



fc2)^k(T)=0. 



(30) 



If fe^ 7^ 0, then behaviour of (fikir) at small times is de- 
termined by the exponent /X3. Obviously, all solutions of 
equation ( pO[) have in general an oscillating character at 
T >~ r,(k) ^ \kz\^/'^'^^~^\ There is no oscillations at 
r ^ T»(k). Among linearly independent solutions there 
exists one, which tends to zero at r — > as ip^'^ ~ r. 

Let us note that the corresponding p'^^ tends to some 
finite value as r — > 0. Another solution, linearly inde- 
pendent on first one, is finite in the zero, but its deriva- 
tive diverges at r ^ as dip"^^ /cLt ~ r^^^'^^, and this 

(2) 

fact means unbounded growth of pj^ and indicates tran- 
sition of the system to the highly nonlinear dynamics 
regime described by the Hamiltonian (|2^). 

For analysis of slow vortical flows, let us write explicitly 
the Hamiltonian 7i*{i7}. For simplicity we put po = 1- 
In accordance with the definition (14) and the expression 
(Eq), we obtain after simple calculations that 



mm = — 

OTT 



1 /•/•dridr2Ea^'^"f^"(ri)^^"(r2) 



Ea^'^° (a:? 



1/2 



(31) 



Let us express the frozen-in field of the vorticity through 
the shape R(i^, ^,t) of vortex hues, (see [D, @, for 
details) : 



(32) 



where v is a. vortex line label, which belongs to the 2D 
manifold M , ^ is a longitudinal parameter along the line. 
The dynamics of vortex lines is determined by the vari- 
ational principle 5Syi — S J Cjidr = 0, with the La- 
grangian of vortex lines having the following form [P, 
1-0: 

CR^^J^d^iy(^(^[R,xR]-R^y^-n,{n{-R}}, 

(33) 

where the subscripts t and ^ mean the corresponding 
partial derivatives. Performing the deformating substi- 
tution 

R= iX,Y,Z) = {r-f^'X, T-f'^Y, r-f'^Z), (34) 
we obtain the Lagrangian for R(i^, ^, r): 



R^ X R 



Rj dfd^i^ 



'Stt, 



(Ra('^i^Ci)-R;2('^2,6)) ^^ ,2 ,2 

— = = atid vidhd V2 

|R(i^l,6)-R(^2,e2)| 



— / [x^YZ{^L2-^Ji:i) 



Y^zx{fi3 - 



Z^XYifii- ^i2)]d^d'ij 



(35) 



It is interesting to note that the same Lagrangian for vor- 
tex lines corresponds to the usual Eulerian incompress- 
ible hydrodynamics in isotropic space with coordinates 
f — (i, y, z) — {t^^x, t 
external straining potential velocity field 



^^y, T^'-^z), but in the presence of 



\p.ix, ^22/, ^^?,z). 



(36) 



In conclusion, let us say a few words about hydrody- 
namics in the homogeneous cosmological models of the 
other types Q, where the space is not flat. In this case, 
a quantitative analysis is difficult due to curvature of 
space and a complicated dependence of metrics on time, 
but the qualitative result about initial free motion regime 
for fluid elements remains valid. The reason is that an 
appropriately deflned relative density p is conserved in or- 
der of magnitude, as well as covariant components Pa of 
the momentum field, while contravariant components 
tend to infinity at i — > 0, like in the flat model. This fact 
means an infinite growth of the quantity PaP°' and ap- 
proaching of the velocity magnitude \/vaV°' to the speed 
of light, regardless field p. Therefore, a weak initial spa- 
tial inhomogeneity in distribution of Pa{9) may result, 
after some time, in a strong spatial inhomogeneity of p. 
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